On the energy landscape of 3D spin Hamiltonians with topological order 
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We explore feasibility of a quantum self-correcting memory based on 3D spin Hamiltonians with 
topological quantum order in which thermal diffusion of topological defects is suppressed by macro- 
scopic energy barriers. To this end we characterize the energy landscape of stabilizer code Hamil- 
tonians with local bounded-strength interactions which have a topologically ordered ground state 
but do not have string-like logical operators. We prove that any sequence of local errors mapping 
a ground state of such Hamiltonian to an orthogonal ground state must cross an energy barrier 
growing at least as a logarithm of the lattice size. Our bound on the energy barrier is shown to be 
tight up to a constant factor for one particular 3D spin Hamiltonian. 
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Topologically ordered phases of matter display a va- 
riety of fascinating properties having no counterpart in 
the classical physics. Most notable ones are topological 
invariants such as the the Hall conductance ground 
state degeneracy 0, and topological entanglement en- 
tropy 0,3 which are insensitive to generic local pertur- 
bations [5|-l7|. The intrinsic stability against perturba- 
tions motivated interest to topological phases as a stor- 
age medium for a reliable quantum memory [&, 9J and 
as a platform for building a topological quantum com- 
puter 

A big open question in the theory of topological quan- 
tum order (TQO) concerns feasibility of a non-volatile, 
or, self-correcting, quantum memory [6|,lll|. Such a mem- 
ory would permit reliable long-term storage of quantum 
information in a presence of sufficiently weak thermal 
noise without need for active stabilization and error cor- 
rection during the storage period. The main challenge 
in designing Hamiltonians with self-correcting properties 
is to combine TQO with an energy landscape that could 
prevent errors caused by thermal fluctuations from ac- 
cumulating. This could guarantee that the error density 
remains sufficiently small during the entire storage period 
and the encoded information can be safely extracted from 
the memory by performing an active error correction at 
the read-out phase. 

In spite of being intrinsically stable against perturba- 
tions at the zero-temperature, TQO models display ex- 
treme fragility against thermal fluctuations 12] suggest- 
ing impossibility of quantum self-correction. A thermal 
stability analysis involving finite-temperature extensions 
of the topological entanglement entropy has been under- 
taken for the 2D and 3D toric code models by Castelnovo 
and Chamon [l^] , and by Iblisdir et al [3] ■ These models 
were shown to undergo a transition from a topologically 
ordered phase at T = to a different phase with either 
partial or no topological order at any positive tempera- 



for the toric code models was performed by Alicki et 
al 16|, |l7|. It showed that the 4D toric code Hamilto- 



181, Chesi et al 



121 

20| 



ture 



The first rigorous analysis of self-correcting properties 



nian has self-correcting properties for sufficiently small 
temperature, while 2D and 3D toric codes are not self- 
correcting at any finite temperature. The ideas of [If 
were developed further by Kay 
and Pastawski et al 12 ll. 



The main feature of the 4D toric code model responsi- 
ble for self-correction is a macroscopic energy barrier th&t 
must be crossed by any sequence of local errors whose 
combined action on encoded states cannot be corrected 
at the final read-out phase Jv^ . The height of this barrier 
grows linearly with the lattice size due to a finite string- 
tension characterizing boundaries of membranes associ- 
ated with errors. It is analogous to the energy barrier 
separating ground states with positive and negative mag- 
netization in the ferromagnetic 2D Ising model. Unfortu- 
nately, this behavior cannot be reproduced in any known 
2D or 3D model due to a presence of point-like excita- 
tions carrying a non-trivial topological charge, or, point- 
like defects 22]. These defects are analogous to domain- 
walls in the ID Ising model — a single isolated defect has 
only a constant energy cost, but its creation requires a 
highly non-local operation affecting a macroscopic num- 
ber of qubits (spins). Whether or not the presence of 
point-like defects rules out self-correcting properties may 
depend on how fast these defects can diffuse across the 
system. For example, Hamma et al used a coupling 
with a bosonic field to create an effective long-range at- 
tractive interaction between defects whereby suppressing 
the diffusion. A different possibility is realized in the 3D 
Chamon's model [13, [l^ . This model offers a topologi- 
cal protection against diffusion of some types of defects 
(but not all of them). These defects, called monopoles 
in [25], can be created at corners of rectangular shaped 
membranes. A hopping of a single isolated monopole 
between adjacent lattice sites is a highly non-local oper- 
ation affecting a macroscopic number of qubits, see [25l | 
for details. 



2 



In the present paper we propose yet another possibihty 
to suppress the diffusion of defects that can be reahzed in 
certain 3D spin Hamiltonians with strictly local bounded- 
strength interactions. These Hamiltonians, associated 
with stabilizer error correcting codes [2^, have a pecu- 
liar property that isolated defects cannot move further 
than a certain constant distance away without creating 
other defects. For brevity, we shall refer to this property 
(stated more formally below) as a no-strings rule because 
it is closely connected to the absence of logical string-like 
operators capable of moving the defects. Let us point 
out that the first example of a 3D spin Hamiltonian with 
TQO obeying the no-strings rule has been found only 
quite recently by one of us (i^. We prove that any se- 
quence of local errors creating an isolated defect from the 
vacuum with no other defects within distance R must 
cross an energy barrier at least c log R for some constant 
c. Furthermore, the length of such error sequence must 
grow at least as R^ for some constant 7 > 1. The same 
bound applies to creation of any isolated cluster of defects 
with a non-trivial total topological charge. 

It shows that although defects do not interact directly, 
their diffusive motion is suppressed by logarithmic energy 
barriers preventing the defects from spreading (the con- 
cept of a diffusive motion must be used with care in our 
case because individual defects can only move a constant 
distance away). 

We also prove a logarithmic lower bound on the energy 
barrier for implementing any logical operator. More pre- 
cisely, we prove that any sequence of local errors mapping 
a ground state to an orthogonal ground state must cross 
an energy barrier at least clogL, where L is the lattice 
size and c is some constant. For the Hamiltonian discov- 
ered in [i^] this bound is tight up to a constant factor. 
Although the scaling of the energy barrier is not as favor- 
able as the one in the 4D toric code, we point out that 
the energy barrier does not grow with the lattice size at 
all for all previously studied TQO Hamiltonians in the 
2D and 3D geometry. A naive estimate of the storage 
time r for a memory with an energy barrier B operat- 
ing at a temperature T can be made using the Arrhenius 
law, namely, r e^^'^. Since in our case B = clogL 
for some constant c, we arrive at r L'^/'^ . Although 
this 'derivation' gives only polynomial scaling of r, the 
degree of the polynomial can be made arbitrarily large 
by choosing sufficiently small temperature. 

It is worth pointing out that 2D Hamiltonians with 
TQO always have string-like logical operators and thus 
3D is the smallest spatial dimension for constructing 
Hamiltonians obeying the no-strings rule. Indeed, it was 
shown by Terhal and one of us [2^ that for any 2D local 
stabilizer-type Hamiltonians the energy barrier for im- 
plementing at least one logical operator is constant. It 
should also be noted that 3D translation-invariant sta- 
bilizer Hamiltonians with TQO can obey the no-strings 
rule only if the ground state degeneracy is not invariant 



under changing lattice dimensions 2^ 3^ 



An obvious difficulty in proving a lower bound on the 
energy barrier is that there are too many ways to choose 
an error path (i.e., a sequence of local errors) implement- 
ing a fixed Pauli operator. In fact, we do not impose any 
restrictions on the length of the error path except that 
it is finite. Our bound should be applicable to any such 
error path. An additional difficulty is that a Pauli opera- 
tor that an error path needs to implement is only defined 
modulo stabilizers. We resolve these difficulties using a 
novel technique which can be regarded as a renormaliza- 
tion group in the space of error paths, see Fig. [1] for a 
brief summary of the technique. Let us now state our 
main results more formally. 




FIG. 1. Renormalization group technique used to prove a 
logarithmic lower bound on the energy barrier. Horizon- 
tal axis represents time. Vertical axis represents RG level 
p = 0, 1, . . . ,Pmax- A sequence of level-0 errors (single-qubit 
Pauli operators) implementing a logical operator P defines 
a level-0 syndrome history (yellow circles) that consists of 
sparse (S) and dense (D) syndromes. The history begins and 
ends with the vacuum (0). For any level p > 1 we define a 
level-p syndrome history by retaining only dense syndromes 
at the lower level. A syndrome is called dense at level p if 
it cannot be partitioned into clusters of size < (10a)'' sep- 
arated by distance > (lOa)''^^, where a is a constant coef- 
ficient from the no-strings rule. Each level-p error (horizon- 
tal arrows) connecting syndromes S' ,S" is equivalent to the 
product of all level- (p — 1) errors between S' , S" modulo a 
stabilizer. We prove that these stabilizers can be chosen such 
that level-p errors act on 2"'-^' qubits. Since at the highest 
level p — Pmax a single level-p error is a logical operator, one 
must have Pmax = i^(log L). We prove that level-p dense syn- 
dromes contain n(p) defects which implies that at least one 
syndrome at level p — Pmax — 2 consists of Q{logL) defects. 



Stabilizer code Hamiltonians. Quantum spin 
Hamiltonians based on stabilizer codes provide conve- 
nient mathematical models of TQO that became an ac- 
tive research topic over last years. We consider a regular 
D-dimensional cubic lattice A with periodic boundary 
conditions and linear size L, that is, A = Each site 
u e A is populated by a finite number of qubits. A large 
class of TQO models can be described by the so-called 
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stabilizer Hamiltonians defined as 

M 



(1) 



a=l 



where each term Ga is a multi-qubit Pauh operator (a 
tensor product of /, X, Y, Z with an overaU ±1 sign) and 
different terms commute with each other. Hence 



GaGh — GfjGa and G^ 



(2) 



The abehan group Q generated by Gi , . . . , Gm is called 
a stabilizer group of the code. Elements of Q are called 
stabilizers. Physically realistic Hamiltonians may only 
involve interactions between small subsets of qubits lo- 
cated close to each other. We shall assume that each 
generator Ga acts non-trivially (by X, Y or Z) only on a 
set of qubits located at vertices of an elementary cube. It 
is allowed to have more than one generator per cube. Any 
short-range stabilizer Hamiltonian can be written in this 
form by performing a coarse-graining of the lattice. The 
Hamiltonian may or may not be translation-invariant. 

We shall assume that H is frustration- free 3l|, that 
is, ground states ipo of H obey Ga ipo = tpo for all a, 
or, equivalently, the stabilizer group Q docs not con- 
tain — /. Consider any multi-qubit Pauli operator E. 
A state tp ~ Eipo is a,n excited eigenstate of H. Obvi- 
ously, Ga ip = iV' where the sign depends on whether Ga 
commutes (plus) or anticommutes (minus) with E. Any 
flipped generator {Ga V' = ""0) will be referred to as a 
defect. Let us emphasize that a configuration of defects 
in ijj is the same for all ground states V'o- An eigenstate 
with m defects has energy 2m above the ground state. 
For brevity, we shall use the term vacuum for a ground 
state of H whenever its choice is not important. A Pauli 
operator E whose action on the vacuum creates no de- 
fects is either a stabilizer {E Cz G), or a logical operator 
{E ^ G, but E commutes with G)- In the former case 
any ground state of H is invariant under E. In the latter 
case E maps some ground state of H to an orthogonal 
ground state. 

Topological order. A Hamiltonian is said to have 
topological order if it has a degenerate ground state 
and different ground states are locally indistinguishable. 
We shall need a slightly stronger version of this condi- 
tion that involves properties of both ground and excited 
states. These properties depend on a length scale Ltqo 
that must be bounded as Ltqo ^ for some constant 
/3 > 0. (For stabilizer codes, Ltqo corresponds to the code 
distance.) Most of stabilizer code Hamiltonians with 
topological order satisfy our conditions with Ltqo ^ L. 
Our first TQO condition concerns ground states: 

// a Pauli operator E creates no defects when applied 
to the vacuum and its support can be enclosed by a 
cube of linear size Ltqo then E is a stabilizer, E E G ■ 



Our second TQO condition concerns excited states. A 
cluster of defects S will be called neutral if it can be 
created from the vacuum by a Pauli operator E whose 
support can be enclosed by a cube of linear size Ltqo 
without creating any other defects. Otherwise we say 
that 5 is a charged cluster. Given a region A C A we shall 
use a notation Br (A) for the r-neighborhood of A, that 
is, a set of all points that have distance at most r from 
A. Here and below we use /oo-distance on Z£. We shall 
need the following condition saying that neutral clusters 
of defects can be created from the vacuum locally (32i] . 

Let S be a neutral cluster of defects and Cmin{S) be 
the smallest cube that encloses S . Then S can be cre- 
ated from the vacuum by a Pauli operator supported 

on Bi{Cynin{S)). 

No-strings rule. We can now state the property of 
having no logical string-like operators. Informally, it says 
that applying an operator with a 'string-like' support to 
the vacuum cannot create charged defects at the end- 
points of the string, assuming that the string is suffi- 
ciently long. Let us now define this property rigorously. 
Let E be any Pauli operator whose support can be en- 
closed by a cube of linear size Ltqo and 5 be a cluster 
of defects obtained by applying E to the vacuum. Let 
Ai , A2 be any pair of disjoint cubes of the same linear 
size p. We shall say that _B is a logical string segment 
with anchor regions Ai , A2 iS S is contained in the union 
AiU A2. Equivalently, E must commute with all genera- 
tors Ga located outside AiU^2- We will say that a logical 
string segment E has aspect ratio a iff the distance be- 
tween Ai and A2 is at least ap. We shall only consider 
string segments with sufficiently large aspect ratio, say, 
a > 1 . A logical string segment E is called trivial iff the 
cluster of defects contained inside any anchor region is 
neutral. 

A code has no logical string-like operators iff there 
exists a constant a such that all logical string seg- 
ments with aspect ratio greater than a are trivial. 

In the rest of the paper we shall abbreviate the condition 
of having no logical string-like operators as a no-strings 
rule. We note that a 3D stabilizer code (Code 1) discov- 
ered in [i^l obeys our topological order conditions with 
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tqo 



L and obeys the no-strings rule with a = 15. 



Energy barrier. Let us consider a process of building 
a logical operator P from local errors. It can be described 
by an error path — a finite sequence of local Pauli errors 
El, ... , Et such that P = Et ■ ■ ■ E2E1. For simplicity 
we shall assume that each local error Et is a single-qubit 
Pauh operator X, Y, or Z. Applying this sequence of 
errors to a ground state ipo generates a sequence of states 
{i^{t)}t=o,...,T, where tp{0) = ipo and tp{T) = PipQ are 
ground states of H, while the intermediate states ip{t) — 
Ef ■ ■ Ei'ipQ are typically excited. We say that a logical 
operator P has energy barrier lo iff for any error path 
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implementing P at least one of the intermediate states 
4>{t) has more than ui defects. Note that we do not impose 
any restrictions on the length of the path T (as long as 
it is finite). In particular, one and the same error may 
be repeated in the error path several times at different 
time steps. We shall also consider an energy barrier for 
creating a cluster of defects S from the vacuum. We will 
say that S has energy barrier w iff for any Pauli operator 
E that creates S from the vacuum and for any error path 
implementing E at least one of the intermediate states 
has more than uj defects. 

Our main results are the following theorems. Both 
theorems apply to any stabilizer Hamiltonian Eq. ([ij on 
a I?-dimensional lattice that obeys the topological order 
condition and the no-strings rule. 

Theorem 1. The energy barrier for any logical operator 
is at least clogL, where L is the lattice size, and c is a 
constant coefficient. 

Theorem 2. Let S be a neutral cluster of defects con- 
taining a charged cluster 5' C S* o/ diameter r such that 
there are no other defects within distance R from S' . If 
r + R< Ltqo, then the energy barrier for creating S from 
the vacuum is at least clogR, where c is a constant co- 
efficient. 

The constant c depends only on the spatial dimension D, 
the constant a in the no-strings rule, and the constant (3 
in the bound Ltqo > ■ Since the proof of both theorems 
uses the same technique, below we shall focus on proving 
Theorem [TJ Proof of Theorem [5] will require only minor 
modifications that are explained later on. 

Proof of Theorem [II A configuration of defects created 
by applying a Pauli operator E to the vacuum will be 
called a syndrome caused by E. The process of build- 
ing up a logical operator P by a sequence of local er- 
rors El, ... , Et can be described by a syndrome history 
{S{t)}t=o.,...,T- Here S{t) is the syndrome caused by the 
product Ef ■ ■ El, that is, the partial implementation of 
P up to a step t. The syndrome history starts and ends 
with the vacuum, i.e., S{0) = S{T) = 0. Without loss 
of generality all intermediate syndromes S{t) are non- 
empty. For any integer p > define a level-p unit of 
length as (33| 

e(p) = (io«)f, p = o,i,.... 

Let S{t) be any non-empty syndrome. Recall that each 
defect in S{t) can be associated with some elementary 
cube of the lattice. 

Definition 1. A syndrome S{t) is called sparse at level 
p iff the set of elementary cubes occupied by S(t) can 
be partitioned into a disjoint union of clusters such that 
each cluster has diameter at most ^(p) and any pair of 
distinct clusters combined together has diameter larger 
than ^{p -\- 1). Otherwise S{t) is called dense at level p. 



For example, suppose all defects in S{t) occupy the 
same elementary cube. Since an elementary cube has 
diameter 1, such a syndrome S{t) is sparse at any level 
p > 0. If S{t) occupies a pair of adjacent cubes, S{t) 
is sparse at any level p > I, however it is dense at level 
p — 0. We note that the partition of S{t) into clusters 
required for level-p sparsity is unique whenever it exists. 

Lemma 1. Suppose a non-empty syndrome S{t) is dense 
at all levels q = 0, . . . ,p. Then S{t) contains at least p+2 
defects. 

Proof. Let C^"^ , . . . , Cg^-* be elementary cubes occupied 
by S{t). Obviously, S{t) contains at least g defects. Since 
S{t) is non-empty and dense at level 0, we have g > 2 and 
there exists a pair of cubes ci°'' , C^°'' such that the union 
ci"'' U C^"'' has diameter at most ^(1). Combining the 
pair into a single cluster we obtain a partition 

S{t) = ^ U . . . U C^i\ where each cluster ci^^ has 
diameter at most ^(1). Suppose S{t) is dense at level 1. 
Then g > 3 and there exists a pair of clusters Ca^"* , C^^"* 
such that the union Ca^^UC^^^ has diameter at most ^(2). 

Combining the pair da\cl^^ into a single cluster and 
proceeding in the same way we arrive at g > p -\- 2. □ 

Let us define a level-p syndrome history as a subse- 
quence of the original syndrome history {S{t)}t=o,...,T 
that includes only those syndromes S{t) that are dense 
at all levels q = 0, . . . ,p — 1, see Fig. [T] The level-0 
syndrome history includes all syndromes S{t). The syn- 
drome history starts and ends with the vacuum (empty 
syndrome) at any level p. Let S{t') and S{t") be a con- 
secutive pair of level-p syndromes. We define a level-p 
error E connecting S{t') and S{t") as the product of all 
single-qubit errors Ej that occurred between S{t') and 
S{t"). Level-p errors are represented by by horizontal 
arrows on Fig. [TJ We would like to show that E can 
be regarded as an approximately local error on a coarse- 
grained lattice characterized by the unit of length ^{p). 
The problem however is that we do not have any bound 
on the number of single-qubit errors Ej in the interval 
between S{t') and S{t"). In the worst case, E could act 
non-trivially on every qubit in the system. The following 
lemma shows that level-p errors can be 'localized' by mul- 
tiplying them with stabilizers. Let m be the maximum 
number of defects in the syndrome history, such that any 
S{t) contains at most m defects. 

Lemma 2. Let S' = S(t') and S" = S{t") be a consec- 
utive pair of syndromes in the level-p syndrome history. 
Let E be the product of all errors Ej that occurred be- 
tween S' and S" . // 4m(2 + S,{p)) < Ltqo, then there 
exists an error E .supported on B^(^p-^{S' U S") .such that 
EE is a stabilizer. 
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The proof of the lemma, presented in Appendix A, 
uses induction in the level p. The proof relies crucially 
on the no-strings rule. The latter asserts that an isolated 
charged defect belonging to some sparse syndrome cannot 
be moved further than distance a away by a sequence of 
local errors. Since the no-strings rule is scale invariant, 
it can also be applied to a coarse-grained lattice to show 
that isolated charged clusters belonging to some sparse 
level-p syndrome cannot be moved further than distance 
a£,{p) away, see Appendix A for details. 

Let Pmax be the highest RG level, that is, the smallest 
integer p such that a single level-p error E maps the vac- 
uum to itself, see Fig.[T] We claim that Pmax = ^i}ogL). 
Indeed, suppose that 'im{^{pmax) -I- 2) < Ltqo- Then 
we can apply Lemma [5] to the level-pmax syndrome his- 
tory with S' = S"' = (vacuum). Lemma [2] would 
imply E = I, that is, E must be a stabilizer. On the 
other hand, E is equivalent to a logical operator mod- 
ulo stabilizers. Hence we obtain a contradiction unless 
im{^{pmax) + 2) > Ltqo- Wc cau assume that the max- 
imum number of defects is m ^ log L (if not, there is 
nothing to prove). Since Ltqo grows as a power of L, we 
conclude that Pmax = ^(^ogL). The syndrome history 
must contain at least one syndrome S{t) which is dense 
at all levels q = 0, . . . ,Pmax — 2 since otherwise Pmax 
could not be the highest RG level. Lemma [T] then im- 
plies that such syndrome S{t) contains i}{\ogL) defects 
proving Theorem [H □ 

We shall prove Theorem [2] using exactly the same ap- 
proach as above. Let S' be a neutral cluster of defects and 
i? be a Pauli operator creating S from the vacuum, with 
S' C S of diameter r being charged. Consider a hierarchy 
of syndrome histories similar to the one shown on Fig. [TJ 
where we now maintain the initial syndrome and the fi- 
nal syndrome S for all levels. Let Pmax be the highest RG 
level. Then a single level-pmax error E creates S from the 
vacuum. Suppose 4TO(^(pmax) -t- 2) < Ltqo, where m is 
the maximum number of defects in the syndrome history. 
Lemma [2] implies that E is the equivalent modulo stabi- 
lizers to E supported on (5*). If ^{pmax) < i?/4, 
then E must act on two separated regions, one near S' 
and another far from S' . This means that S' alone can 
be created by a Pauli operator whose support is enclosed 
by a cube of linear size r -|- R/2. Since r + R < Ltqo, it 
is contradictory to the assumption that S' is charged. 
Therefore, £,ipmax) > ^/4, or p^ax = ^{logR). In 
case where 4:m{^{pmax) + 2) > Ltqo > Li, we also have 
Pmax = ^^(logi?) provided m -C logi?. (There is nothing 
to prove if m ~ logi?.) Since there must be a syndrome 
that is dense for all level-p where p — 0,1, . . . ,Pmax — 2, 
Theorem follows from Lemma [T] 

Optimality. Our lower bound on the energy barrier 
is optimal up to a constant factor. As we noted earlier. 
Code 1 in 27[ exhibits TQO with Ltqo ^ L and obeys the 
no-strings rule with a constant a = 15. This code can be 



IZ ZI IX XI 

/ \ / 

XI- 



ZI — -^zz 



II- 



■IZ 



II 



XX 



■IX 



IZ ZI IX XI • o 

FIG. 2. Stabilizer generator for Code 1 in [131 . This 
translation-invariant model exhibits TQO, and obeys no- 
strings rule with a — 15. The diagram on the right is a 
cube in the dual lattice. The filled dots indicate the defects 
created by XI at the center of the cube, which we call level-0 
pyramid. 



described by a stabilizer Hamiltonian Eq. ([T]) with two 

qubits per site and two stabilizer generators per cube 

as shown on Fig. [5] We shall use a notation PQu for 

a two-qubit Pauli operator P ® Q applied to the pair 

of qubits located at site u. A single-qubit error XI^ 

creates a cluster of 4 defects (violated generators of Z- 

type) located at a set of cubes s'f '' = {c, c+x, c+y, c+z}, 

where the center of c is obtained from u by a translation 

— {x + y+z)/2, see Fig. [2] We shall refer to S'c°^ as a level- 

pyramid with an apex c. For any integer p > define 

a level-p pyramid with an apex c as a cluster of 4 defects 

located at a set of cubes S^^'' = {c, c+2Px, c+2Py, c+2Pz}. 

We note that a level- (p-|- 1) pyramid with an apex c can 

be represented as a sum (modulo two) of four level-p 

pyramids with an apex at c, c-\-2Px, c + 2Py, and 0-1-2^2. 
(p) 

Therefore, Sc can be created from the vacuum by an 
error Eu defined recursively as 

~ ^u+2Px^u+2Py-^u+2Pz 

with Eu'^ = XIu- Simple induction shows that E'if^ acts 
on 4^ qubits. Its support can be regarded as a fractal 
object with a fractal dimension 2. 

Suppose that the lattice has periodic boundary con- 
ditions and its linear size is L = 2" for some integer n. 
Then the four defects of a level-n pyramid cancel each 
other, i.e., ^c"^ = 0. It shows that the error Eu ^ is 
either a stabilizer or a logical operator. Consider an aux- 
iliary operator Z„ = Yii j=i Zlu-x+iy+jz- One can easily 
check that commutes with all stabilizer generators of 
X-type, see Fig. [51 On the other hand, Z„ anti-commutes 
with since their supports overlap only on a single 
site u -\- (2" — l)x = u — X and their restrictions on this 
site {ZI and XI respectively) anticommute. It shows 
that i?!"-* is a logical operator. Let us now show that 

(n) 

E'u' has energy barrier at most 4 -t- 41ogL. 

We claim that the error E^'' creating a level-p pyramid 

(v) 

Sc has energy barrier ujp — 4p + 4. Below we shall only 
consider error paths implementing single-qubit factors of 

E^'^ in a certain order. We shall use induction to prove 

(p) 

the following statement: Eu can be implemented by a 
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sequence of single- qubit errors without creating more than 
ojp defects. The apex c contains a defect after each error 
in this sequence. For p = 0, the claim is obvious. Sup- 
pose we have aheady proved the claim for p = 0, . . . , g. 
Let us build a \evel-{q + 1) pyramid Sc'^'^^^ . We start 
from building a level-g pyramid Sc'^\ which requires at 
most 4(7 + 4 defects and the apex c contains a defect af- 
ter each elementary error. Now we can build another 
level-g pyramid S'^'^2ix using at most 4g -I- 4 defects plus 
the defects that were already present in Sc'^\ Although 
si'^^ has 4 defects, one of them is located at the apex of 
^c'+2ix- "^^^ induction hypothesis implies that any par- 
tial implementation of s''^^2ix ^^^'^ contains a defect at 
c -I- These two defects cancel each other and must 
be subtracted from the total cost. Thus we can build two 
nearby pyramids at cost of 4 -I- (4g + A) — 2 = Aq-\- 6 de- 
fects. Both pyramids together have 6 defects. The same 
argument shows that building the third pyramid S^^^2iy 
requires at most 6 + {Aq -|- 4) — 2 = 4(7 -|- 8 defects, and the 
three pyramids together have 6 defects. Building the last 
pyramid S^^^^iz requires at most 6 + {Aq + 4) — 2 = 4(7 + 8 

defects. Thus, we have constructed a pyramid s'f'^^'^ us- 
ing at most 4(g -I- 1) -I- 4 defects. Note that the apex 
c was occupied by a defect over all stages of this con- 
struction. It proves the step of induction. We conclude 

in) 

that a logical operator Eu has energy barrier at most 
4n + 4 = 4 -I- 41ogL. 

Separation of charged clusters. We have shown 
that in a process of isolating a charged cluster, there 
is a logarithmic energy barrier. The following theorem 
quantifies how long the process must be. 

Theorem 3. Let E he a Pauli operator creating S , a neu- 
tral cluster of defects containing a charged cluster S' Q S 
of diameter r such that there are no other defects within 
distance R from S". If r + 2R < Ltqo, then the weight of 
E must be > cR'' for some constant 7 > 1 and c. 

The proof again makes use of renormalization group, 
and shows that there is a subset of 'fractal dimension' 
7 > 1 in the support of E. We assume that E has weight 
minimum possible. Let be a odd positive number. We 
say a set of sites C C A is a level-p chunk if diam(C) < 
wP . A path in the lattice is a finite sequence of sites 
M2, • ■ • , Un) such that d{ui,Ui+i) = 1. (Recall that 
we use the ^oo metric d.) Equipped with paths, we can 
say whether a set is connected. 

Definition 2. A connected level-p chunk C Q S is 
maximal with respect to a set of sites S if there ex- 
ist a connected subset C'° ^ C and a path C 
(ui, . . . , m, . . . , u„) C C° satisfying 

(i) d{ui,Un) ^wP- 

(ii) d{ui,m),d{un,m) > ^'''^Jf'' — ^ 



(Hi) C° contains the connected component of m in 

S„p_„p-i (to) n S, and 
2 

(iv) C contains the connected component of C° in 

B^P-i {c°)ns. 

2 

The last two conditions restricts the position of C in C 
such that lies sufficiently far from the boundary of C. 
The site m will be refereed to as a midpoint of C. Let S 
be the support of the Pauli operator E, any restriction 
of which obeys no-strings rule. 

Lemma 3. Given a path in S joining ui and Un such 
that d{ui,u„) = lujP — 1, there are I disjoint maximal 
chunks of level-p whose midpoints are on ^. 

Lemma 4. For sufficiently large w, a maximal level- 
{p-\-\) chunk C with respect to S admits a decomposition 
into w -h 1 or more maximal chunks of level-p with respect 
to S. 

The proofs of Lemma [3l|4] can be found in Appendix B. 
The support of the minimal Pauli operator E in Theo- 
rem [3] must admit a path connecting S' and S\S' . Oth- 
erwise, S' can be regarded as being created locally, and 
our topological order condition demands the cluster be 
neutral. Since the path has length > R, Lemma [3] says 
we have a maximal chunk of level-p where p is such that 
ujP < R < w^"*"^. Lemma[3]implies any maximal chunk of 
level-p must contain at least {w -I- 1)^' sites. This proves 
Theorem [3] with 7 = '"fp™^"'""' > 1 . A similar argument 
proves the lower bound d = ^l{L'^) on the code distance d 
of Code 1 in since the minimal logical operator must 
contain a path of length L. 



CONCLUSIONS AND OPEN PROBLEMS 

We proposed an argument in favor of self-correcting 
properties for a class of 3D spin Hamiltonians with topo- 
logical quantum order. Thermal diffusion of defects in 
these Hamiltonians is suppressed by the presence of log- 
arithmic energy barriers. A novel technique for proving 
lower bounds on the energy barriers is presented. 

Our results rise several questions. Firstly, Theorem [2] 
implies that the energy landscape contains a macroscopic 
number of local minimums separated by macroscopic en- 
ergy barriers. These minimums correspond to low-lying 
excited states with only a few defects such that separa- 
tion between defects is of order L. Such energy landscape 
suggests a possibility of a spin glass phase at sufficiently 
low temperature. We note that a spin glass phase can 
indeed be realized for some classical spin Hamiltonians 
with logarithmic energy barriers such as the model dis- 
covered by Newman and Moore 



34| . Interplay between 



topological order and spin glas siness has been studied re- 
cently by several authors 2J, [35[ . Secondly, our paper 
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leaves open the question of how to perform error cor- 
rection necessary to extract encoded information from a 
memory. In particular, no efhcient error correction algo- 



rithm is known for the stabilizer code discovered in [27 1 . 
One could speculate however that locality and a macro- 
scopic distance of the code are sufficient conditions for 
having a non-zero error threshold under random indepen- 
dent errors. Finally, an exciting open question is whether 
any 3D stabilizer code Hamiltonian with TQO has point- 
like topological defects. A proof of this conjecture would 
give an evidence that "strong self-correction" similar to 
the one in the 4D toric code is impossible for realistic 
physical systems. 
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APPENDIX A 

In this section we prove Lemma [2] For convenience of 
the reader, we repeat the statement of the lemma. 

Lemma 2. Let S' = S{t') and S" = S{t") be a consec- 
utive pair of syndromes in the level-p syndrome history. 
Let E he the product of all errors Ej that occurred be- 
tween S' and S" . If 4m(2 + ^(p)) < Ltqo, then there 
exists an error E supported on i3^(p)(5" U S") such that 
EE is a stabilizer. 

Proof. Let us use induction in p. The base of induction 
is p = 0. In this case E = Ej is a single-qubit error. If 
the qubit acted on by E does not belong to Bi{S' U S"), 
one must have S' = S". It means that E is a single- 
qubit error with a trivial syndrome. Topological order 
condition implies that E is a stabilizer. Choosing E = I 
proves the lemma for p = 0. 

Suppose we have already proved the statement for 
some levelp and let us prove for levelp-|-l. Let S' = S{t') 
and S" = S{t") be consecutive syndromes in the level- 
(p -I- 1) history. Consider first the trivial case when 
S' — S{t') and S" — S{t") are also consecutive syn- 
dromes in the level-p history. Then 5" and S" are con- 
nected by a single level-p error E which, by induction 
hypothesis, has support on Sj(p) (S" U S"') (modulo stabi- 



lizers). The latter is contained in B^i^p^i-f {S' U S") which 
proves the induction step. 

The non-trivial case is when there is at least one level-p 
syndrome between S' and S" . The interval of the level-p 
syndrome history between S" and S" can be represented 
(after properly redefining the time variable t) as 



S' 



^5(1) ^5(2) 



£2, 



Here all syndromes S{1), . . . , S{t) are sparse at the level 
p and all transitions are caused by single level-p errors. 
The condition of sparsity implies that the set of elemen- 
tary cubes occupied by S(t) has a unique partition into 
a disjoint union of clusters Ca{t) such that each cluster 
has diameter at most ^(p) and the distance between any 
pair of clusters (if any) is at least 

dist(a(i),Cfc(i)) >e(p+l)-2e(p) > (10«-2)e(p) 
> 8<(p). 

Represent any intermediate syndrome as a disjoint union 



S'(i) =5"(t)U5"(0, t = l,.. 



(3) 



where S'^{t) and S'"(t) include all charged and all neu- 
tral clusters Ca{t) respectively. Let g be the number of 
clusters in S'^{t). We claim that g does not depend on t. 
Indeed, since a level-p error Et acts on ^(p)-neighborhood 
of S{t) U S{t -\- 1), the sparsity condition implies that Et 
cannot create/annihilate a charged cluster Ca{t) from the 
vacuum, or map a charged cluster to a neutral cluster (or 
vice verse). The same argument shows that each cluster 
Ca{t) C 5^(i) can 'move' at most by ^(p) per time step, 
that is, we can parameterize 

S^{t) = Cl{t)U...UCg{t) 

such that a 'world-line' of the a-th charged cluster obeys 
the continuity condition 



dist{Cait+l),Ca{t)) <aP)- 



(4) 



We can now use the no-strings rule to show that all 
charged clusters are 'locked' near their initial positions, 
so that their world-lines are essentially parallel to the 
time axis. More precisely, we claim that 

dist(Ca(t), Ca(l)) < a^ip) for aU 1 < i < r. (5) 

Indeed, suppose Eq. ([5]) is false for some a. Using the 
continuity Eq. ([4]) one can find a time step ti such that 
dist(Ca(ti),Ca(l)) > aiip) and dist(Ca(t), C,(l)) < 
Q^(p) for all 1 < t < ti. Let Edose be the product 
of all level-p errors Ej that occurred between S{1) and 
S{ti) within distance (2 -I- a)^(p) from Ca(l). Since all 
intermediate syndromes are sparse at level p, the net ef- 
fect of Eciose is to annihilate the charged cluster Cq(1) 
and create the charged cluster Ca(ii). Equivalently, ap- 
plying Eciose to the vacuum creates a pair of charged 
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clusters Ca(l) and Ca(<i). However, this contradicts to 
the no-strings rule since Ca(l) and Ca{ti) have linear size 
at most ^{p) while the distance between them is greater 
than a£,{p). Thus we have proved Eq. ([S]). 

We will say that af S A is close to S' iff x G ;B^(p+i) (5'). 
We will say that a? e A is close to S" iff f G B^(p+i)(S"'). 

Let Et be the level-p error causing the transition from 
S{t) to S(t + 1), where t = 1, . . . , r - 1. Let is the 
restriction of Et onto Sj(p)(5=(i) U S%t + 1)), and E^ is 
the restriction of Et onto U + 1)). The 

sparsity of S{t) and localization of level-p errors then 
implies that 

Et = E^ ■ El\ (6) 

We claim that any error E^ is close to S' . Indeed, each 
cluster in 5'^(1) is within distance 2^{p) from S' since 
otherwise a single level-p error Eiead would be able to 
create a charged cluster from the vacuum. Using Eq. ([5]) 
we infer that Ca(t) G B(^2+a)^{p)iS') for all a = 1, . . . , 
Therefore E^ is close to 5". 

We shall now define a 'localized' leading error Eiead 
that maps the syndrome S' to 5^(1) such that the sup- 
port of Eiead is close to S'. For any neutral cluster 
C G S'"(l) let 0'(C) be a Pauli operator creating C from 
the vacuum. By definition of a neutral cluster, we can 
assume that 0'{C) has support in Bi{C). Set 

Eiead ~ Eiead ■ G ■ 0'(C), 

ces"(i) 

where G G ^ is a stabilizer to be chosen later. Since 
Eiead is a single level-p error, we have two options: 

(a) any cluster C G is within distance ^{p) from 
5', 

(b) there is exactly one cluster Cfar G such that 
the distance between Cfar and 5" is greater than 

In case (a) we set G = / if the support of Eiead is close 
to S' and G = E'iead otherwise (if Eiead is not close to S", 
but any cluster G G S'"(l) is within distance £,{p) from 
S", the error Eiead creates no defects when applied to the 
vacuum, that is, Eiead must be a stabilizer). In case (b) 
the error Eiead creates a neutral cluster Cfar from the 
vacuum. In this case we will set G = EieadO' [C far)- We 
conclude that the support of Eiead is close to S' and Eiead 
maps S' to S"=(l). 

We can apply the same rules to the error Etaii and the 
syndrome S'"(t). It allows us to find a stabilizer G and 
define a localized error 

Etaii — Etaii ■ G ■ Y[ 0"{G), 

such that Etaii maps 5^(t) to 5" and the support of 
Etaii is close to S". The operator 0"(G) above creates a 
neutral cluster G G 5'"(t) from the vacuum. 



We can now define a localized level- (p -|- 1) error E 
whose support is close to S' U S" as 

E = Eiead ■ El - ■ ■ E'^_i - Etaii- (7) 

By construction, it describes an error path 

S' S\\) ^ S%2) • . • S^t) S". 

It remains to check that E - E is a stabilizer. Combining 
Eq. ([6]) and Eq. ^ we conclude that 

E - E ^ {EieadElead) ■ E" - - - • {EtailEtail) - 

Applying E - E to the vacuum generates the following 
chain of transitions: 

vac ^'-°''^'-°''> 5"(1) ^"(2) 

■ - - 5"(t) ^'"'^'°"> vac. (8) 

Here all transitions are caused by errors whose support 
can be enclosed by at most m cubes of linear size 2 + ^(p). 
Each syndrome S"{t) consists of at most m neutral clus- 
ters of diameter ^{p), i.e., it can be created from the 
vacuum by an error whose support can be enclosed by 
at most m cubes of linear size 2 + (,{p)- Now the state- 
ment that _E • £' is a stabilizer follows from the following 
proposition. 

Proposition 1. Let Qj be Pauli operators causing a 
chain of transitions 

vac — > bi — > D2 — > - - - — > br > vac. 

Let Pj he some Pauli operator creating the syndrome Sj 
from the vacuum. Suppose the support of any operator Pj 
and any operator Qj can be enclosed by at most n cubes 
of linear size R such that AnR < Ltqo- Then the product 
Q = Qi - - - QrQr+i is a stabilizer. 

Proof of Proposition [1], Let tpo be any ground state. 
Define a sequence of states 

^(1) = PiQi • V'o, 

i^{j + i)^{p,p,+i)Q,+i-m, 

^{r + l) = Qr+iPr-i^{r), 
for = 0, . . . , r — 1. Obviously, 

V'(j) = Pj • (Qi • • • Q]) • V'o for j = 1, . . . , r 

and ilj(r + 1) — Q - ijjQ. It follows that all states are 
ground states, and the transition from ■(/'(j) to + 1) 
can be caused by a Pauli operator 

Gj = PjPj+iQj+i- 

Let Mj be the support of Gj. Obviously, Mj can be en- 
closed by at most 3n cubes of linear size R. Suppose 
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first that Mj is a connected set, i.e., one can connect any 
pair of qubits from Mj by a path . . . ,ui) such that 
the distance between Ua and Ua+i is 1. Then Alj can 
be enclosed by a single cube of linear size at most 3nR. 
Our assumptions on n and i?, and the topological order 
condition then implies that Oj is a stabilizer. Suppose 
now that Mj consists of several disconnected components 
M " , such that the distance between any pair of compo- 
nents is at least 2. Let O" be the restriction of Oj onto 
a connected component M". Locality of the stabilizer 
generators Ga and the fact that Oj commutes with any 
Ga implies that O" commutes with any Ga- Further- 
more, Af" can be enclosed by a cube of linear size at 
most 3nR. Topological order condition then implies that 
O" is a stabilizer. It follows that all operator Oj are sta- 
bilizers, that is, + 1) = for all j. However, it 
means that Q'tpo — tpo for any ground state V'Oi that is, 
Q is a stabilizer. □ 

□ 

APPENDIX B 

In this section we prove the technical lemmas needed 
for Theorem [3l 

Lemma 3. Given a path in S joining ui and u„ such 
that d{ui,Un) = Iw^ — 1, there are I disjoint maximal 
chunks of level-p whose midpoints are on 

Proof. For convenience, we assume that the z-coordinates 
of ui and u„ are and Iw'' — 1, respectively. Con- 
sider / -|- 1 planes Pi perpendicular to the z-axis, whose 
z-coordinates are iw^ for i = 0,1,..., L In each 
region between the two consecutive planes Pi-i and 
Pi, there is a subpath Q = (uj._j^, . . . ,Uj.) such that 
^("ji-ii '^ji) ~ — 1. Choose mi G Q such that 
d{uj._-^,mi), d{mi,Uj.) > ^''^'^ ■ Let Gi be the connected 
component of within S'nSmP (m^). Add, if necessary, 

some points of Q to Ci to get a maximally connected G^. 
This C- is a maximal chunk of sites with midpoint being 
nii. Any two Cj"s are disjoint since each of them lies in 
a unique region enclosed by Pi-i and Pi. □ 

Lemma 4. For sufficiently large w, a maximal level- 
(p+l) chunk G with respect to S admits a decomposition 
into w + 1 or more maximal chunks of level-p with respect 
to S. 

Proof. Recall that S is the support of the Pauli operator 
E, any restriction of which obeys no-strings rule. Define 
the boundary of a subset U oi S to be dU = Bi (U) HUT] 
S. Then, any subset U of sites with boundary enclosed in 
a two disjoint regions can be regarded as a string segment. 

By the definition of the maximal chunk, there ex- 
ists a path (wi, . . . , m, . . . , w„) in G° C G such that 



d{ui, Un) = wP~^^—wP. We assume that the z-coordinates 
of ui, Un differ by w^"*"^ —w^. We will show that there are 
sufhciently long and separated paths in G° , to which we 
apply Lemma [3] to find w + 1 maximal chunks of level-p. 
They will lie in ^?™£.(C°), and hence in G. 

Let M (N) be the subset of S consisted of sites whose 
z-coordinates differ from that of ui by at most rjut^. 
First, suppose dG° is not contained in M U TV. Since 
ui E M and m„ S N, there is a site s & G° adjacent 
(of distance 1) to dG° such that d{s, ui), d{s, Un) > rjw^ . 
Furthermore, o?(s, to) > "'^ 2~™'' ~ Otherwise, G° con- 
tains a site in the boundary, which is a contradiction. 

Consider the shortest network M of paths in G° con- 
necting four sites ui,m,Un, s. (The length of a network 
of paths is the number of sites in the union of the paths.) 
Let C be the shortest path in Af from ui to u„. If s is 
not contained in B^^piO^ then C' ^ A/" joining s to a site 
on C has a subpath ^" C (' of diameter at least 2wP such 
that C" is separated from (hy vu^. Applying Lemma [3] to 
C and C", we find at least w + I maximal chunk of level- 
p. If TO is not contained in B^^piOi ^ similar argument 
reveals at least w + 1 maximal chunk of level-p. 

Suppose both s and to are contained in BswpiO- Ob- 
serve that C \ {B4wp(s) U B4u]p{m)) consists of three con- 
nected components Cii C27 Csi two of which have diameter 
> Sw^ and the other has diameter > (ry— 4)u)P. 

Two distinct Bu^iCi) and B^iQ) {ij = 1,2,3) do 
not overlap because of the minimality of C. Applying 
Lemma [31 we find w + rj ~ 21 maximal chunks of level-p. 
Choosing 77 > 21, we get the desired result. 

Next, suppose dG° is contained in M U N. Let 
SM,SN e G° \ {M U N) be sites adjacent to AI and 
N, respectively. The separation between M and N is 
{w — 1 — 2r])w^. If it is greater than rj' aw^ , there must 
be a z-plane P that contains sm or sjv such that P n C° 
has diameter > rj'w^; Otherwise, the no-strings rule is 
violated. Let vi,V2 G P n C° be sites separated by rj'w^. 
The four sites, wi, m„, vi,V2 are sufficiently separated and 
connected by some paths in G° . Arguing as before, we 
find (w— 1)+77' — 16 maximal chunks of level-p. The choice 
of rj' > 17 and w > 1 + 2ri + r/'a proves the lemma. □ 
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